Quasi-Periodic Oscillations in Black Hole Candidates as an Indicator of 
Transition Between Low and High States 

Robert E. Rutledge"'^'^, Walter H. G. Lewin^, Michiel van der Klis^'"^, Jan van Paradijs'^'^, Tadayasu 
Dotani^, Brian Vaughan^, Tomaso Belloni'^, Tim Oosterbroek^, and Chryssa Kouveliotou* 

^ Department of Astronomy, University of California, Berkeley, CA 94720 

^ 37-627, Dept. of Physics, Massachusetts Institute of Technology, Cambridge, MA 02139 

^ Astronomical Institute "Anton Pannekoek", University of Amsterdam, Center for High Energy Astrophysics, Kruis- 
laan 403, 1098 SJ Amsterdam, The Netherlands 

* Physics Department, University of Alabama in Huntsville, Huntsville AL 35899 

^ Institute of Space and Astronautical Science, 3-1-1 Yoshinodai, Sagamihara Kanagawa 229-8510, Japan 
® Space Radiation Laboratory, MS 220-47, California Institute of Technology, Pasadena, CA 91125 

Astrophysics Division, Space Science Department of ESA, ESTEC, P.O. Box 299, 2200 AG Noordwijk, The Nether- 
lands 

* Universities Space Research Association, Huntsville, AL 35800 



ABSTRACT 

By comparing positions on a spectral color-color diagram from 10 black hole candidates 
(BHCs) observed with Gmga (1354-64, 1826-24, 1630-47, LMC X-1, LMC X-3, GS 2000-f25, 
GS 2023-f33, GS 1124-68, Gyg X-1, and GX 339-4) with the observed broad-band noise (0.001- 
64 Hz; BBN) and quasi-periodic oscillation (QPO) variability, we find that the so-called "Very 
High State" is spectrally intermediate to the Soft/High-State and Hard/Low-State. We find 
a transition point in spectral hardness where the dependence of the BHC QPO centroid fre- 
quency (of GS 1124-68 and GX 339-4) on spectral hardness switches from a correlation to an 
anti-correlation; where the BBN variability switches from High-State to Low-State; and where 
the spectral hardness of the QPO relative to that of the BBN variability is a maximum. This 
coincidence of changing behavior in both the QPO and the broad-band variability leads us to 
hypothesize that the QPO is due to interaction between the physical components which dominate 
the behaviors of BHCs when they occupy the Hard/Low and Soft/High States. We conclude that 
these QPO should be observed from BHCs during transition between these two states. Compar- 
ison with QPO and BBN behavior observed during the 1996 transition of Cyg X-1 supports this 
hypothesis. 

We also report 1-3 Hz QPO observed in GS 2000-1-25 and Cyg X-1 in the Hard/Low State, 
and we compare these to the QPO observed in GS 1124-68 and GX 339-4. 



Subject headings: black hole physics - stars: oscillations -X-rays: stars ~ stars: Cyg X-1 - stars: 
GX 339-4 - stars: GS 1124-68 



1. Introduction 



The X-ray (2-20 keV) phenomenology of galactic black hole candidates (BHCs) has been widely 
studied (see Tanaka & Lewin 1995 and Van Der Klis 1995 for recent reviews). The BHCs have many 
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similarities in their phenomenology, in that similar X-ray spectral hardnesses are accompanied by similar 
fast-timing behavior (0.01-100 Hz). Often, it is this similarity of behavior alone which identifies an individual 
X-ray source as a BHC, when this source has no measured mass-function, and it presents correlated X-ray 
spectral and timing behavior which is similar to that observed from BHCs which do have a (high, ^^3 Mq) 
measured mass-function. 

This X-ray phenomenology is largely composed of three "states" , historically tied to the observed 
source intensity, which has evolved to be tied to the source photon spectral hardness. The "Soft/High" state 
(HS; historically, when the 1-10 keV source flux was high) is marked by a soft photon spectrum dominated 
by an "ultra-soft" black-body component, with kT $1 kev, and weak (~ few %) rms (root-mean-square) 
variability which is a power-law as a function of frequency, of slope ~ 1. The "Hard/Low" state (LS) is 
marked by a hard photon spectrum dominated by a power-law component, and constant %rms below some 
frequency (^ 0.01-10 Hz), decreasing roughly as a power-law above this frequency. The "Very High" state 
(VHS) is marked by the presence of quasi-periodic oscillations (QPO) with frequencies between 3-10 Hz. In 
addition, an "Off" state (in which the source is quiescent) is generally included in any phenomenology. More 



recently, an "Intermediate State" (IS; Belloni et al. 1996; Belloni et al. 1997; Mendez fc Van Per Klis 1997) 



has been recognized, in which the variability is band-limited, constant below some break frequency 0.1 — 10 
Hz) and decreasing above this, while the energy spectrum shows both an ultra-soft black-body component 
and a power-law component >10 keV. The distinction made between this state and the HS and LS are 
the simultaneous presence of the ultra-soft and power-law energy spectral components, with band-limited 
broad-band noise (BBN) variability and (in GS 1124-68) the presence of QPO at source fluxes below that of 
the so-called "High State". 

As a structure for phenomenologically characterizing BHCs, the state paradigm lacks quantitative 
value. For example, Cyg X-1 may be observed in the Soft/High State. One might well ask, how high? Since 
"high" describes the observed flux, the uncertainty in distances to BHCs makes it difficult to inter-compare 
the behavior of objects on this basis. There currently exists no gradation of states, no formal quantitative test 
to determine if a particular source occupies a state or not, or "how much" it occupies a state. The description 
of additional "states" , without somehow tying them through some invariant observational parameter to other 
"states" does not illuminate the relationship between the various states. In addition, it seems reasonable 
that the BHCs do not occupy "states" at all, but that BHC behavior changes slowly as a function of its 
observables, since individual BHCs are observed to have a wide range of luminosities, spectral hardnesses 
and power density spectra, rather than three or four of these. As such, investigating BHC behavior as a 
function of continual variables may be a useful observational approach. 

In an individual source, observational states are defined on the basis of qualitative changes in its 
power-spectrum and energy-spectrum, which are observed to be strongly correlated. It seems likely that the 
instantaneous accretion rate (M) is the underlying variable parameter responsible for the observed changes 
in source behavior. However, between different sources, other parameters - such as the mass of the black 
hole, or the stellar type of its companion - may also come into play in determining the observed source 
intensity behavior. Comparisons in behaviors between different objects therefore may eventually provide 
observational handles on these parameters. 

The combined use of color-color diagrams (CCDs) and the contemporaneous timing behavior - 
usually indicated by the strength of source variability and its functional dependence on frequency in Power 
Density Spectra (PDS) - has quite successfully revealed regularities in the complex QPO and BBN timing 
behavior of GS 1124-68 ( [Takizawa et al. 1997| ) and GX 339-4 ([Miyamoto et al. 199lD . Following this work, it 



was found that in the VHS, in which the QPO are found, the PDS (apart from the QPO) can be like those in 
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the Hard/Low State, or like those in the Soft/High ( Miyamoto et al. 1993| ; Takizawa et al. 1997). In a similar 
vein, it has been shown that, as the fraction of the total intensity due to the power-law spectral component 



increases above 10%, the %rms variability increases, as well ( Miyamoto et al. 1994 ); the variability behavior 
of GS 1124-68 was demonstrated to be correlated to the strength of its hard-spectral-component, which had 
been long suspected. 

The CCD vs. fast-timing s approach was first used in X-ray astronomy to untangle correlated timing 



and spectral/intensity behavior of non-pulsing neutron-star Low-Mass X-ray Binaries (LMXBs; [Hasinger fc 



Van Per Klis 1989 ). Consequently, it has become clearer what the accretion rate is when the various distinct 
behaviors of these objects are observed. More detailed theory on the spectral and timing behavior has been 
produced, as well as proposals for the relationship between the behavior of neutron-star LMXBs and that of 



BHCs and pulsars (Van Der Klis 1994a ,b ; see Van Der Klis 1995 for review) 



We have completed a study of the simultaneous spectral and timing behavior of ten BHCs observed 
with Gmga (1354-64, 1826-24, 1630-47, LMC X-1, LMC X-3, GS 2000-f 25, GS 2023-1-33, GS 1124-68, Cyg X- 
1, and GX 339-4). A description of BBN behavior of these sources is being presented elsewhere (Rutledge 



et al. 199'^ ; hereafter BHCLFN). Here, we describe and discuss the spectral and energy dependence of QPO 
in the four sources which exhibit them in the Ginga data (GS 1124-68, GX 339-4, Cyg X-1 and GS 2000-f 25) 
and the strength and frequency dependence of the underlying BBN. These results are a summary of analyses 



presented in greater detail elsewhere ( Rutledge 1996 ) 



We have examined the source spectral state, as characterized by its location on a CCD, and simulta- 
neous timing properties. We find that BHCs show QPO typically when they have a spectral hardness ratio 
between the Hard/Low State and Soft/High State. When we follow the timing behavior of the two objects 
which exhibit this QPO (GX 339-4 and GS 1124-68) as a function of spectral hardness, changes in behavior 
of the QPO and the underlying BBN indicate that the so-called VHS is a transition between the Hard/Low 
State and the Soft/High State. This has important implications for the origin of the QPO, as an indicator 
of this transition, perhaps being produced by interaction between two distinct physical processes responsible 
for the HS and LS phenomena. 

An important note on terminology: Historically, "Low State" was applied to describe the behavior 
of BHCs when the source intensity was low, and "High State" for when it was high. We are faced with a 
descriptive problem, in that the behavior of the BHCs we describe herein - aside from intensity, and the 
presence of QPO - is indistinguishable from that observed from the "Low State" , or the "High State" . It is 
a common aspect of theories of BHC accretion to ascribe the Low State to one physical component, such as 
a geometrically thick accretion disk or a spherical corona, and the High State behavior to another, such as 



a geometrically thin accretion disk {e.g. 


Narayan & Yi 1994; 


Abramowicz et al. 1995; 


Narayan & Yi 1995|; 


Chen et al. 1995 


; Chen & Taam 1996; 


Narayan 1996; 


Esin et al. 1997). Therefore, throughout this paper we 



use the terms "Low State" and "High State" not to indicate source intensity, but behavior which is ascribed 
to the respective physical components we mention above. 



2. The X-ray Data and Analyses 



All X-ray data on black hole candidates (BHCs) in the present work were obtained with the Large 
Area Counter (LAC) on the the Ginga satellite ( Turner et al. 1989 ). The LAC on-board Ginga had a 
geometric area of ^ 4000 cm^, a nominal energy range of 1-37 keV, divided into 48 energy channels, with 
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energy resolution of 18% (i?/6keV)^/^. The time resolution of the data, depending on the data-acquisition 
mode, was typically between 0.001-16 sec. 



2.1. Color-Color Diagram 



The data used in the CCD are described in detail elsewhere (BHCLFN; App. B in Rutledge 1996 ). 
Briefly, our data are drawn from all data acquired by Ginga during observations of each of ten X-ray black 
hole candidates (1354-64, 1826-24, 1630-47, LMC X-1, LMC X-3, GS 2000-^25, GS 2023-f33, GS 1124-68, 
Cyg X-1, and GX 339-4), which had energy resolution sufficient to produce our hardness ratios (see also 
below). All data were corrected for background, deadtime, and aspect. Following these corrections, we 
produced hardness ratios using three bands - the Low Energy Range (2.3-4.6 keV), the Medium Energy 
Range (4.6-9.2 keV) and the High Energy Range (9.2-18.4 keV). These bands were combined into the Hard 
Ratio (High Energy Range / Medium Energy Range) and the Wide Ratio (High Energy Range / Low Energy 
Range), which we use to produce the CCD. Each point represents 128 seconds of contiguous integration time. 
Only points which have error bars in both colors showing significantly measured hardness ratios (>4cr) are 
shown (the remainder are discarded). Due to inadequate background subtraction, not all data acquired by 
Ginga on a particular source is used in this CCD. Thus, the CCD is biased against fainter intensities, for 
which either the background subtraction techniques employed or counting statistics in a 128-sec integration 
were inadequate. After all selections were applied, a total of 680 ksec (of 1500 ksec) of data from all sources 



were used in the CCD. The CCD is presented in Sec. 3.1 



In the lowest energy bin, the effect of galactic absorption is minimized (although not eliminated) by 
using a high lower energy bound (2.3 keV). All BHCs in the present study have column densities between 



1.4 10^1 cm-^and 6.0 lO^i cm-^ except GS 2000+25 (8.4 lO^^ cm^^) and 1630-47 (25 lO^^ cm^^; Van 



Paradijs 1995 ). The effect of galactic absorption in the 2.3-4.6 keV energy range is small - for a black-body 
of kT=0.5 keV, the Wide Ratio changes by <10% between 1.4 lO^^ cm'^and 6.0 lO^i cm'^; at 25 lO^^ 
cm~^, this becomes a 30% effect in the 2.3-4.6 keV energy range. 

Throughout the present work, we refer to spectral "states" , typically as indicated by the contem- 
poraneous Wide Ratio value. For brevity, we use a progression of letters (A, B, C, D..., J) which indicate 
a range of Wide Ratio space, each half a decade wide in logarithmic space, increasing in sequence with in- 
creasingly hard spectrum (A is —4.0 <Wide Ratio < —3.5, B is —3.5 <Wide Ratio < —3.0, etc.). The width 
of these ranges was arbitrarily chosen, with the idea that the source behavior does not vary significantly 
within a single range. This is not always the case, particularly when QPO are observed, and it then becomes 
necessary to further subdivide the data based on the spectral values into "sub-states" {e.g. Fl, F2, ...). 
However, the letters A-J, are retained for convenience in referencing the timing analysis to a simultaneous 
spectral indicator. 



2.2. Timing Analysis 

In Table |^, we show the quantity of data which are available to us for the timing study, separated 
by source and time resolution. In all, 853 ksec are available for the fast-timing analysis; of this, 305 ksec 
have time resolution of or better than 0.06250 sec. 
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An extensive description of our PDS production, including mathematical definitions, normaliza- 



tion, and computation of correction factors for dead time and background, is to be found elsewhere (Rut- 



ledge 1996). We repeat the pertinent details here. 



2.2.1. Composite PDS 

To make efficient use of all available data, each observation was used to produce any (or all) of three 
possible types of FFTs, which we define here: 

• Low Frequency FFT: Data is binned to 16.0 s, with a total duration of 1024 s, using 64 bins of data. 
The resulting frequency range is 0.0009765-0.03125 Hz. 

• Medium Frequency FFT: Data is binned to 0.0625 s, with a total duration of 64 s, using 1024 bins 
of data. The resulting frequency range is 0.015625-8 Hz. 

• High Frequency FFT: Data is binned to 7.8125 msec, with a total duration of 8 s, using 1024 bins 
of data. The resulting frequency range is 0.125-64 Hz. 

Thus, a single observation with a duration of 1500 s and a time resolution of 7.8 msec is used to 
produce 187 High-Frequency FFTs, 23 Medium-Frequency FFTs and 1 Low-Frequency FFT. 

Using raw data selected from those listed in Table |l] (listed in greater detail in App. B of Rutledge 
1996 ), we produced FFTs in three separate energy bands: the Low Energy Range (=2.3-4.6 keV), Medium 
Energy Range (=4.6-9.2 keV) , High Energy Range (=9.2-18.4 keV) and in the summed Intensity Range 
(=2.3-18.4 keV). This produced a library of Low-, Medium-, and High-Frequency FFTs in three independent 
photon energy ranges and one summed photon energy range from which we draw to produce PDS, based on 
the selection criteria we apply for any given analysis. 

To produce a Composite PDS, individual FFTs are selected by the Wide Ratio, measured during 
the time period covered by the FFT. Each individual FFT is made into a single PDS, its calculated dead- 
time corrected Poisson Level subtracted, and the resulting power corrected for background and instrumental 
deadtime. Finally, all the resulting PDS are averaged together into a single PDS. The uncertainty in the 
power at each frequency is taken to be the measured standard deviation in the power divided by the square 
root of the number of measurements. Thus, the uncertainty in the power contains any (non-biased) systematic 
uncertainty in the Poisson-subtraction, dead time and background %rms correction factors. 

These PDS are logarithmically rebinned, and overlapping frequency ranges are dropped in the PDS 
with the lesser integrated time. For example: with 50 Medium-Frequency FFTs and 200 High Frequency 
FFTs, the Medium-Frequency FFTs represent 64 s x 50 = 3200 s of data, while the High-Frequency FFTs 
represent 8 s x 200 = 1600 s of data. Therefore, all the data <8 Hz would be dropped in the High-Frequency 
PDS, and then the two PDS are merged. 

The final frequency range of the composite PDS is 0.00097-64 Hz or less if data of sufficient con- 
tiguous duration or time resolution are not available, and is normalized as rms'^/Hz, and corrected for dead 
time and background. 

The data of the high, medium and low frequency range parts of the single PDS are not necessarily 
taken at the same time. Thus, this type of PDS is subject to systematic uncertainties due to changes in the 
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PDS within a single selection criteria. In addition, there may be more data available in one frequency range 
than in others. 

For each of the ten individual objects, we produced four composite PDS (of the Low, Medium, and 
High energy Range, and for the full Intensity Range) from data selected by Wide Ratio to be in each of the 
coarse spectral hardness bins (A, B, C, ...), as illustrated in Fig. |l| and described in Sec. 3J. We used these 



PDS to investigate the dependence of BBN variability on spectral hardness, to search for the presence of 
QPO and understand the dependence of QPO properties on spectral hardness and photon energy. 

There are, at times, intensity "spikes" which occur in the data, a few msec in duration, which are 
typically not due to source variability, but instead to a variety of background sources (some of these spikes 
were due to geomagnetically trapped particles from orbiting nuclear power generators. These events, and 
the observed 2-37 keV PHA spectra, are discussed further elsewhere ( Rutledge 1996| )). We systematically 



excluded these spikes from our analysis, as they can produce significant, measurable power across all frequen- 
cies. When a single time bin is found in which the number of counts is 20 a (calculated from the intrinsic 
scatter of the data used for that particular FFT) greater than the mean number of counts, that value was 
replaced by a value randomly selected from a set of bins in the data used for that FFT. Such events, even if 
real, are identifiable as discrete events, which are not our primary focus. 



2.3. Parameterized Fits to the PDS 

In practice, the models adopted to fit a PDS may not be ideal, or even formally acceptable according 
to statistical tests. Discrepancies may be systematic - for example, due to the skew of a quasi-Lorentzian peak 
which is not accounted for in the model - or they may be due to excursions from simply parameterizable 
continuum behavior in the PDS ~ that is, deviations from the model which are non-biased, like random 
scatter, but statistically significant. There exists no a 'priori reason that such biased and noisy PDS should 
not exist. Thus, while it is desirable to fit the models as closely to the data as possible, there can be a trade- 
off between the accuracy and simplicity of a model. These parameterizations are not physically motivated, 
and thus their utility is to accurately describe PDS using a small set of numbers. 

For the present study, we desire simple parameterizations which adequately describe PDS from 
all 10 sources. We avoid the introduction of new parameters, apart from switching to a more reasonable 



parameterization. We use a maximum-likelihood fitting method (least Press et al. 1995), and quote the 
Icr single parameter uncertainties, the reduced x^, and the number of degrees of freedom. 

We use a number of models, most of which have become widely used in describing PDS of X-ray 
binaries, to fit the observed power (P) as a function of frequency [y) to the PDS. The model shapes we used 
to fit to the real source variability are all continuous functions of frequency, and include: 

• Power-Law, Pc>c v~'^^ 

• Flat-Top Single-Power-Law, a constant power below some frequency v\ breaking to a Power-Law with 
slope a\ above this frequency. 

• Broken Power-Law, represented by one power-law slope (ai) below a critical frequency V\ and a 
different, typically steeper, power-law with slope a2- If we find a\ ^0.3 (that is, if the first power-law 
is very flat), then we call this model "Flat-Top Single-Power-Law" instead of a "Broken Power-Law". 
This is a definition of convenience, having no theoretical or observational impetus. 
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• Flat-Top Broken- Power-Law, a constant power below some frequency vi breaking to a Power-Law with 
slope ai above this frequency, breaking again to a (typically steeper) Power-Law a2 above frequency 

1^2, with Vi < V2. 

• Double Flat- Top Single-Power-Law, a model in which the power is constant below some frequency 
i^i, breaking to a power-law of slope ai, which evens out into a constant value again at some higher 
frequency i^2(> t^i), and then breaks again into a power-law of the same slope (ai) at some higher 
frequency v^. 

• Flat-Top Power-Law Exponential, a model in which the power is constant below some frequency vi, 
and continuous with the product of a power-law and an exponential function above this frequency. 

In fitting models to the data, we employ one of the above models, which we call the continuum 
model, or Broad-Band Noise (BBN) model, and a model to account for quasi-periodic oscillations from: 

• Lorentzian, used to model quasi-periodic oscillations, which has a centroid frequency Uc and Full- Width 
at Half-Maximum (FWHM; F), and a normalization related to the %rms variability. 

• Harmonics, sometimes added to account for the appearance of Harmonics and sub-Harmonics to the 
QPO; the frequencies and FWHM values of these Harmonics are constrained to be 0.25, 0.5, 2.0 and 
4.0 that of the Fundamental, except where noted. 

The %rms variability of the Lorentzians and Harmonics are integrated between —oo and oo Hz. 

2.4. Energy Dependency of the PDS 

We present the results of the fits to the four sources in which we observed QPO (GS 1124-68, GX 339- 
4, Cyg X-1, and GS 2000+25) in four energy ranges. For each PDS, we indicate the mean source spectrum 
by the Wide Ratio (=High Energy Range/Medium Energy Range). Finally, we fit the PDS with BBN and 
QPO models to determine the parameterization of the PDS and the %rms variability in each spectral state. 

3. Results 

3.1. Color-Color Diagram 

In Fig. 0, we show a color-color diagram (CCD) of all X-ray data in the present study, showing the 
Wide Ratio vs. the Hard Ratio. 

The dashed-line divisions on the figure are artificial boundaries, separating the spectral parameter 
space into equally sized areas based on the Wide Ratio. These divisions are labelled with the alphabetic 
shorthand (A, B, C, D, etc.), encompassing two such ranges, corresponding to a decade in magnitude of the 
Wide Ratio. 

The approximate locations of each source are labelled on the CCD. The CCD shows that, for most 
of these sources, spectral evolution in these colors produces roughly parallel tracts which run diagonally 
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across the CCD. These tracks show smooth variation in the spectrum within the CCD, although there are 
discontinuities ~ hkely due in part to short observing times relative to the time-scale of spectral evolution. 

We have approximately indicated the positions on the CCD where BBN behavior expected from the 
various "states" was observed in the PDS of sources in that area. Sources which are spectrally soft, have 
PDS which have a single-power-law dependence on frequency, and present only weak (~few %rms) variability 
(1-20 keV) occupy the area labelled "High State". Sources which are spectrally hard, have complex PDS - 
typically, but not always, constant below a frequency '^0.1-10 Hz, and decreasing as a power-law at higher 
frequencies - and have strong ('^20-50%rms) variability (1-20 keV) occupy the area labelled "Low State". 
The behavior of High State and Low State variability will be presented in detail elsewhere (BHCLFN). 
Sources which present QPO occupy the spectral range labelled "Very High State" . 

Some sources were only observed to occupy the Low State (such as GS 2023-1-33, 1630-47, 1826-24, 
and Cyg X-1). Others were observed only in the High State (1354-64 and LMC X-3). Some were observed 
to traverse two or more of these states (GS 1124-68, GX 339-4, and LMC X-1 - although the observation of 
LMC X-1 in the spectral region corresponding to the Very High State did not have sufficient time resolution 
to confirm the presence of QPO). 

In Fig. ^, we show on the same scale as Fig. |l| the points at which QPO were identified (by eye) as 
present in the PDS of several sources in the present study. In addition, we include a point from QPO from 



GX 339-4 in a previous study Grebenev et al. 1991, observed with the ART-P spectrometer, with a centroid 
frequency of 0.8 Hz, FWHM of 0.15 Hz, and a %rms variability of 7.0±0.8% in the 3-25 keV energy range. 
When this QPO was detected from GX 339-4, the energy spectrum was power-law, with a photon slope of 
1.733±0.004. This places GX 339-4 at (Wide Ratio, Hard Ratio) vector of ( 0.32, 0.40), found by propogating 
the observed photon spectral model through the Ginga response (see Sec. ^.l.l| ). While the Wide Ratio of 
the sources in the present study span four orders of magnitude, the QPO is largely clustered within <1 order 
of magnitude, located close to a single locus in the CCD. Sources which were observed to exhibit QPO were 
also observed at other times at different positions on the CCD, but with no QPO. 

Fig. H shows the Wide Ratio for each of the ten sources (same data as Fig. |l|, described above). There 
are three sources - GS 2000-1-25, LMC X-1, and LMC X-3 - in which some observed Wide Ratio values are 
similar to those observed from GX 339-4 and GS 1124-68 when 3-10 Hz QPO were observed (the Very High 
State spectral range). Timing analysis of all observations (but one) of GS 2000-1-25 produced upper-Hmits to 
QPO strength well above that observed from GX 339-4 and GS 1124-68; in the one observation near Wide 
Ratio=0.1 (at the edge of the Very High State spectral range), 2.6 Hz QPO from GS 2000-1-25 was observed 
(see Sec. |3.2.4|) . LMC X-1 and LMC X-3 were observed with time-resolution of 16 sec - too long to measure 



QPO at 3-10 Hz. Thus, we do not observe any BHCs which exhibit spectral hardness in the Very High State 
spectral range, but no QPO (with upper limits below the observed strengths in GX 339-4 and GS 1124-68). 



3.1.1. Intrinsic Source Spectrum to Wide Ratio Conversion 

As the Wide Ratio (as measured by Ginga) is the parameter to which we tie the observed behavior, 
we provide a figure to relate this ratio to the observed photon spectral parameters of a power-law, black- 
body, and the relative strength of the two. To do so, we simulated photon energy spectra with black-body 
temperatures of kT =0.1, 0.25, 0.4, 0.55, 0.7, 0.85, 1.0, and 1.4 keV; with photon power-law spectra of 
slopes = 0.9, 1.3, 1.7. 2.1, 2.5, and 2.9. (We use the black-body spectrum as a simple parameterization, 
although more physically reasonable models, such as a disk black-body spectrum, may more accurately 
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describe the spectrum in this energy range; Mitsuda et al. 1984; Makishima ct al. 1986| ). We combined these 



two spectral components with different relative strengths, indicated by the ratio (X) of 2-20 keV photon flux 
(photons/sec) of the power-law to that of the black-body, for values of X in the range of 0.0001-100. 

The relationships between X and Wide Ratio for these assumed photon spectra are shown in Fig. ^. 
To have a Wide Ratio ^0.01, a BHC must be dominated by the black-body (that is, ultra-soft) component. 
To obtain a Wide Ratio -JO.l, BHCs with photon spectral slopes < 2.2 must be dominated by the power-law 
photon spectral component. At steeper spectral slopes, and in between 0.01-0.1 where most of the QPO in 
GS 1124-68 and GX 339-4 are observed, the Wide Ratio is a complex function of X, kT, and the photon 
spectral power-law slope. 



3.2. Power Density Spectra 

Here, we present the parameters describing the QPO and their energy dependence, including the 
strength of the simultaneously observed BB variability and the models used to describe them. The de- 
tailed results pertaining to BBN variability when QPO are not observed, and the detailed parameters and 
energy dependence of the BBN variability when QPO are observed are presented elsewhere (BHCLFN; 
Rutledge 19"9^). 



3.2.1. Cyg X-1 

Of the measured Wide Ratio values for Gyg X-1, 80% were in the range 0.13-0.51, which is within 
the Low-State described in Sec. b.ll 



When Cyg X-1 exhibited QPO, the values of its Wide Ratio are largely in one of two ranges: 
0.12-0.15 and 0.16-0.18. These in turn correspond to two different observational periods (90/129-131 and 
91/156-157, respectively). We sub-divided these data into two spectral groups: Wide Ratio=0.1-0.15 (Gl), 
and 0.15-0.32 (G2). The QPO parameters from the best-fit models are given in Table |2[ 

In Fig. ^, we show the PDS of the Gl and G2 spectral ranges, with the power renormalized as vP. 
For all energy ranges, the PDS are identical above 10 Hz. Below 10 Hz, the PDS are different, and appear 
to be complex functions of frequency. Most evident is the presence of a QPO peak near 1 Hz in the PDS of 
the Gl data. The best fit BBN variability model, without a QPO Lorentzian near ~1 Hz, has a > 13, 
which decreases to ~1.6 when the Lorentzian is added. The %rms variability of this QPO decreases with 
increasing energy range from 11.7±0.5 in the Low Energy Range to 8.2±0.6 in the High Energy Range. 

We also fit a Lorentzian to the G2-data PDS. In the Intensity Range a peak was found at 3.14 Hz, 
with FWHM 0.7±0.3 Hz. These values were fixed for the Low Energy Range, Medium Energy Range, and 
High Energy Range PDS, in which the %rms was consistent with being colorless. 

This QPO peak is very broad (Q ^^1), which would usually result in its being called a "noise" feature, 
as QPO is a term traditionally reserved for higher Q signals. We attempted to distinguish this feature from 



the QPO of GX 339-4 and GS 1124-68 (see Sec. 3^) based on the QPO spectral hardness as a function of 
the source photon spectral hardness, but found we were unable to do so on this basis. We therefore forward 
the hypothesis that this feature is related to the QPO observed in the data from GS 1124-68 and GX 339-4. 
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3.2.2. GX 339-4 

Of the measured Wide Ratio values for GX 339-4, 80% were in the range Wide Ratio=0. 044-0. 56. 

We closely examined the data obtained in the September 1988 observations of GX 339-4, divided 
by source spectrum, which is the only period during which the Wide Ratio was in the F spectral range, and 
during which the VHS QPO was observed. In Fig. ||, we show the Wide Ratio as a function of time during 
these observations. The source spectrum (measured by Wide Ratio) hardens between t=0 and 2 days; over 
the next two days, the source spectrum globally softens, while exhibiting short-term reversals - variously 
hardening and softening. 

We divided the data by spectrum, and labeled the spectral regions (Fl-4; all in the F-state) on the 
hardness-intensity diagram (HID) in the upper-right panel of Fig. |6[ We recognize at least three separate 
tracks in the HID - a nearly constant count-rate track (near 1=5800 c/s 2.3-18.4 keV, crossing the regions 
labeled "F3" and "F4"), a nearly constant Wide Ratio track (near Wide Ratio=0.045 - labeled as "F2"), 
and a steeply sloped track (the Wide Ratio region labeled "Fl"). These tracks correspond to those found in 



Fig. 5 of Miyamoto et al. ( |l99lD and Fig. 5.28 of Ebisawa (|1991|) 



We used the Medium- Frequency and High- Frequency FFTs in each spectral region to produce Com- 
posite PDS, in all energy ranges. We show these PDS in Fig ^, in the bottom 16 panels. The scale of all 
PDS is given in the bottom left panel. When spectrally hardest (F3 and F4), the source shows band- limited 
noise, breaking at a frequency of '^few Hz, in addition to QPO which is strongest in the 9.2-18.4 keV energy 
band. When spectrally softer (Fl and F2), the PDS show less variability, and a freqency dependence which 
is more like a power-law, or perhaps band-limited with a break-frequency near 0.02 Hz. 

We fit the Composite PDS produced from the Fl-4 data of GX 339-4. The models used were Single- 
Power-Law (for Fl), Double-Flat- Top-Single-Power-Law (for F2, with the low-frequency flat-top vi frozen 
to below 0.01 Hz), and Flat-Top-Single-Power-Law (states F3 and F4). For all these, we also fit Harmonic 
QPO, allowing the %rms normalization in the first two harmonics and sub-harmonics to vary. Typically, 
the High Energy Range PDS did not constrain the models well; for these PDS, we froze the BBN model 
parameters (except the %rms) at their best fit values found from the Intensity Range PDS. 

In Table ^, we list the QPO and BBN %rms values found from these fits. In all PDS, the strength 
of QPO increases with increasing photon energy. The strength of the QPO docs not show a monotonic 
dependence on spectral hardness of GX 339-4. Looking at the PDS, one might be misled to believe that the 
%rms of the QPO in the Low Energy Range PDS increases with increasing spectral hardness; our model fits 
do not bear this out. One is apparently misled by the relative strength of the BBN power, which accounts 
for a dominant portion of the %rms at the QPO frequency when the spectrum is hard, but considerably less 
when the spectrum is soft. While we do see dependence between the QPO %rms in the Low Energy Range 
it is not monotonic, nor pronounced, with the %rms values in the range of 1.1-2.9%. 

The QPO FWHM and centroid frequency exhibit some dependency (again, non-monotonic) on 
spectral hardness. In the Intensity Range the FWHM increases from 0.5, to 1.0, to 3.6, and back down to 
2.3 as the spectrum hardens between states Fl-4. The centroid frequency goes from 5.8, to 6.3, to 7.4, to 
6.3 as the spectrum hardens between states Fl-4. In all states, the QPO parameters (centroid frequency, 
FWHM) are energy independent. 
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3.2.3. GS 1124-68 

Of the measured Wide Ratio values for GS 1124-68, 80% were in the range 0.00098-0.38. 

The QPO properties in the summed 1.2-18.4 keV range exhibited by GS 1124-68 have been examined 



in detail (Takizawa et al. 1997, hereafter, T97). We here re-examine these properties in our three energy 



ranges, as well as the Intensity Range, using the Wide Ratio spectral indicator. 

Over short periods of time hours) there is typically very little evolution in the source spectral 
hardness, less than 0.5 dex in hardness ratio (see T97, BHCLFN). We therefore produced PDS over single 
observational epochs. We adopted the same time-periods used by T97, and the same labels for these time 
periods (by number 1-15), with two exceptions: 1) we did not distinguish between observation 5L and 
5H (where "L" stands for periods within the trough of a "flip-flop" - see T97 for a description of this 
phenomenon), while "H" stands for the elevated count-rate periods around the flip-flop) - they are joined 
under the single header "5"; also, we split observation 11 into 11a and lib. The time limits for the data 
used here are listed in Table |4[ along with the frequency range of the resultant PDS, and the mean Wide 
Ratio for each observation period. The uncertainty listed in the Wide Ratio is the la spread of observed 
128-second values (the relatively small values demonstrate that the amount of spectral evolution during a 
single observational epoch is typically small). 



As found previously (T97, and references therein; Belloni et al. 1997 ), QPO is detected during 
observations 1, 2, 3, 5, 10, 11, and 15. The QPO during observation 11 appears is strongest in the second 
half (lib). 

In Table ||, we list the centroid frequency, FWHM, and %rms of the fundamental QPO in each PDS 
in which QPO was identified by T97. In addition, we find (weak) QPO in Observation # 8 and 11a. The 
fundamental was identified for each observation as that QPO peak with the greatest %rms in the 9.2-18.4 
keV range, or, if none are identified, in the 4.6-9.2 keV range. 

The QPO was fit as a series of harmonics, with frequencies and FWHM values which are factors of 
0.25, 0.5, 1.0, 2.0, and 4.0 of those of the fundamental. Typically, only the sub-harmonic at 0.5 i^c and the 
second harmonic at 2 Vc were invoked (permitting the best-fit value of the %rms in each to be determined). 
For a few, however, three or four harmonics have significant %rms values. These harmonics were studied in 



detail elsewhere (T97; Belloni et al. 1997; Rutledge 1996). 



In all cases, the %rms of the fundamental increases with increasing photon energy (except observation 
#15, where the S/N is too low to permit detection in the High Energy Range). 

In addition to the detailed parametric fits to the PDS of the individual observations, we also produced 
PDS using the Fl-4 division done for GX 339-4, for qualitative comparison with that source. These are shown 
in Fig. ^, produced similarly as for GX 339-4 (see Fig. |^). When the source is spectrally hard (F3 and F4), 
the PDS show band-limited noise with a break frequency of a few Hz, along with QPO which is strongest 
in the 9.2-18.4 keV energy band. When the source is spectrally soft (Fl and F2) the PDS shows weaker 
variabilty, and a frequency dependence which is either power-law (F2) or band-limited noise with a break 
frequency of ~0.1 Hz. This behavior is similar to that observed from GX 339-4 using the same spectral 
hardness division. 
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3.24- GS 2000+25 

Of the measured Wide Ratio values for GS 2000+25, 80% were in the range 0.021-0.23. The values 
of the Hard Ratio are typically higher than those observed for other objects with similar Wide Ratio values. 

In a check of PDS of individual observation periods, a period of significant variability was found on 
day 252 of 1988 (Sept 8, ~140 days into the outburst), between 15:45 and 17:24 UT. During this period, 
the source was in a spectral state with Wide Ratio=0. 09-0.1. We show the PDS, from an average of 44 
Medium-Frequency FFTs, in Fig. ||. We fit the data with a Flat-Top-Power-Law function with an additional 
Lorentzian component to account for QPO near 2 Hz. We find that the QPO component is required by 
the data in the Low Energy Range (compare 3.4 for 29 degrees of freedom, or dof, without the QPO 
component to 0.7 for 26 dof with the QPO component). We fit the Flat-Top-Single-Power-Law plus QPO 
model to the Low Energy Range, Medium Energy Range, and High Energy Range PDS made from data 
taken during the period when the LS QPO was observed. The best-fit model parameters are listed in Table ^. 
The centroid frequency (Intensity Range PDS) is 2.63±0.06, FWHM is 1.7±0.2, and the %rms variability is 
consistent with being independent of energy. 

We produced a PDS using data with a Wide Ratio consistent with that observed from GS 1124-68 
and GX 339-4 while they exhibit VHS QPO (F range). Visual inspection of the PDS from this observation 
finds no QPO peak. After finding the best-fit BBN models in the Intensity Range we added a Lorentzian, with 
i'c=4, 5, and separately 6 Hz, with FWHM=1.0 Hz, to obtain 3cr upper-limits of the %rms variability of such 
component, which were 7.2%, 4.8%, and 7.5% (respectively). In the High Energy Range the corresponding 
3cr upper-limits were 29%, 26%, and 30%. As these limits are greater than the measured strength of VHS 
QPO in GX 339-4 and GS 1124-68, we cannot exclude the presence of VHS QPO in this data. 



3.3. The Energy Dependence of the QPO 

In Fig. §, we show the centroid frequency and FWHM of the QPO in GX 339-4 and GS 1124-68 as 
a function of spectral hardness. In both GX 339-4 and GS 1124-68, the centroid frequency is dependent on 
spectral hardness. The centroid frequency is greatest near Wide Ratio=0.055 (±0.01 approximately), and 
decreases as the source becomes spectrally harder and as the source becomes spectrally softer. The same is 
true for the FWHM of the QPO - it is greatest at Wide Ratio=0.055, and smaller as the source is spectrally 
harder or softer. We also include the QPO of Cyg X-1 and GS 2000-f 25 in this figure, which show roughly 
consistent behavior. 

In Fig. |l^, we show the %rms variability as a function of photon energy for QPO and BBN variability 
in GS 1124-68 and GX 339-4. In both GS 1124-68 and GX 339-4 (panels a and d), the VHS QPO %rms 
(averaged over all observations) increases with increasing energy. The BBN average power (panels b and e), 
however, increases from the Low Energy Range to the Medium Energy Range; it then remains constant or 
decreases between the Medium Energy Range and High Energy Range. 

In GS 1124-68, when there is no QPO observed (panel c), the BBN average power increases with 
increasing photon energy across all three energy ranges. To produce this panel, we used observations 6, 7, 
9, 12, and 14, in which VHS QPO were not observed, but the source spectrum was similar to that during 
times when VHS QPO were observed (in the range 0.01-0.03 vs. 0.013-0.155). In the Low Energy Range 
and Medium Energy Range the average broad-band power is considerably larger when QPO is present than 
when QPO is absent (compare 9.8±0.1 and 18.3±0.5 %rms with 2.8±0.05 and 10.0±0.2 %rms). However, 
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in the High Energy Range, the aveage power is comparable (15.4±0.8 %rnis with QPO, 17.1±0.4 %rms 
without). The average ratio of BBN variabihty in the High Energy Range vs the Low Energy Range when 
QPO is absent is 6.1 ±0.2 vs. only 1.6 ±0.1 when QPO is present. This may imply that the presence of 
QPO supresses BBN variability in the high energy range. 

In Fig. we show the energy dependence of the LS QPO in Cyg X-1 and GS 2000+25. The QPO 
in these sources is slightly anti-correlated with energy (in Cyg X-1) or constant (in GS 2000+25). This is 
quite a different energy dependence than observed in the QPO of GX 339-4 and GS 1124-68 (sec Tables ^ 
& H, which are always spectrally hard. In particular, observations #1, 2, and 3 of GS 1124-68, which have 
spectral hardness similar to Cyg X-1 and GS 2000+25 here, are not spectrally flat or soft, but spectrally 
hard. This may imply a different mechanism for these QPO of Cyg X-1 and GS 2000+25 from those of 
GS 1124-68 and GX 339-4. (We are unable to characterize the LS QPO of GX 339-4 reported previously by 
Grebenev ( |1991[ ), as no information on the energy dependence of the variability is available). 



3.4. The Variability Ratio 

To compare the spectrum of the QPO relative to that of the broad-band variability in the sources 
Cyg X-1, GS 2000+25, GX 339-4, and GS 1124-68, we define the variabihty ratio: 

_ %rmsQpo(9.2- 18.4keV) . %rmsBB(9.2 - 18.4keV) 
Rqpo,bb - %rmsQPo(2.3-4.6keV) ^ %rmsBB(2.3 - 4.6keV) ^ ' 



which is a ratio of two ratios, each being the ratio of the High Energy Range %rms variability to the Low 
Energy Range %rms variability, in the QPO and (separately) the broad-band components. Values of the 
variability ratio which are consistent with 1.0 mean the QPO and BBN components have the same energy 
spectrum; values which are greater than 1.0 mean the QPO is spectrally harder than the BBN component. 
The variability ratio is independent of the source spectrum, in the sense that if the spectrum changes while 
the %rms variability as a function of frequency remains the same, the variability ratio also remains the same. 

For GS 2000+25 and Cyg X-1, the spectrum of the QPO is the same as, or slightly softer than the 
spectrum of the BBN variability. The variability ratio of GX 339-4 is significantly higher, which indicates 
that the VHS QPO is spectrally harder than the BBN variability. The same is often true for GS 1124-68: 
the VHS QPO is spectrally harder during observations 1, 3, 5, 10, and lib. During observations 2, 8, and 
15, however, the QPO spectral hardness is within a factor of 2 (spectrally harder) of the spectral hardness 
of the BBN variability. 

In Fig. [ij, we show the variability ratio as a function of the Wide Ratio, for the four objects from 
which we identify QPO. For GS 1124-68 and GX 339-4, the value of the ratio is high (>4) while the Wide 
Ratio is in the range 0.03-0.08. While GX 339-4 was not observed with QPO outside of this range, GS 1124- 
68 was observed with spectral hardness both above and below this range, where the ratio is observed to be 
low (^2). For Cyg X-1 and GS 2000+25, from which QPO were observed in their LS, the value of the ratio 
was observed to be low ($2). The similarity in the value of the ratio observed from GS 1124-68, Cyg X-1, 
and GS 2000+25 while they have similar spectral hardness supports the argument that the QPO of these 
objects are of similar origin. The QPO centroid frequencies of GS 1124-68 were 3-5 Hz (Observations #1-3; 
TablesH and ||), which are low, but not identical to the ~ 1 Hz frequencies of Cyg X-1 and GS 2000+25. 
Fig. O illustrates the correlation between the variability ratio and the QPO centroid frequency measured 
from these sources. 
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Discussion and Conclusions 



4.1. Wide Ratio of BHCs when QPO are Observed 



While there is no simultaneous spectral data for the six hour observation period during which 0.08 



Hz QPO were found in LMC X-1 (Ebisawa et al. 1989), in the 13-hrs following this period (observed with 
16-sec time resolution), the Wide Ratio was in the range 0.033-0.087 (Fig. |) which is very close to the 
Wide Ratio of GS 1124-68 and GX 339-4 when QPO were observed. Although High-Frequency PDS were 
produced by Ebisawa et al. , no 3-10 Hz QPO were reported. However, the QPO can turn on and off on a 



timescale of ~1 sec ( Miyamoto et al. 1991 ; Takizawa et al. 1997 ), correlated with a sudden change in source 
spectral hardness (and intensity - the so-called "flip-flops"), and this may explain the absence of QPO in 
LMC X-1. 

We conclude that we observe no BHCs with spectral hardness 0.02 $ Wide Ratio ;$ 0.1 which do 
not exhibit QPO; consequently, we expect that BHCs with spectral hardnesses in this range should produce 
QPO similar to that described here. 



4.2. The Very High State as a Transition between the Low State and the High State 

Previous investigations using the present data concluded that both band-limited noise variability 
behavior (like that seen in the Low State) and power-law noise variability (like that seen in the High State) 



are observed at various times simultaneous with QPO in GX 339-4 (Miyamoto et al. 1991) and GS 1124-68 
(T97). In both studies, it was shown that these behaviors are observed when the object was relatively 
spectrally hard (in the case of Low State- like noise) and spectrally soft (in the case of High State- like noise). 
Neither investigation commented upon the relationship between this transition and the expected transition 
between the Low and High States. 

In GS 1124-68 (T97), it was also found that the QPO centroid frequency had (at least) two distinct 
functional dependencies on the 2-20 keV intensity, one which was valid when power-law (HS-like PDS) 
variability was observed (in Obs #11 and #15), the other when band-limited (LS-like PDS) variability was 
observed. Comparison with the present results shows that the change in the dependence of QPO frequency 
on intensity observed by T97 occurs at Wide Ratio~0.055 ("the transition point"). 

Because: 



1. the Very High State of GX 339-4 and GS 1124-68 (i.e. when QPO are observed) is spectrally between 
the Low State and High State (0.02 <Wide Ratio <0.2); 

2. power-law PDS, similar to that expected from the spectrally soft High State, are observed when QPO 
are present in both GX 339-4 and GS 1124-68 when these sources are spectrally softer than at the 
transition point (0.02 <Wide Ratio^O.055); in addition, the QPO centroid frequency is correlated with 
spectral hardness in this spectral region; 

3. band-limited PDS, similar to that expected from the spectrally hard High State, are observed when 
QPO are present in both GX 339-4 and GS 1124-68 and when these sources are spectrally harder 
than at the transition point (0.055 $Wide Ratioi:0.2); in addition, the QPO centroid frequency is 
anti-correlated with spectral hardness in this spectral region; 
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it seems natural to interpret the VHS as the transition between accretion dominated by the spherical corona 
(the Low-State) and accretion dominated by the accretion disk (the High-State). The change at the (so- 
defined) transition point in the QPO spectral hardness relative to the BBN variability, as well as the de- 
pendence of the QPO centroid frequency on the Wide Ratio from a correlation to an anti-correlation, may 
indicate that the QPO is a manifestation of the interaction between the corona and the accretion disk. 

It is widely assumed that the accretion rate (M) is correlated with the observed flux in BHCs, 
and that the 2-20 keV flux is used as a monotonic indicator of this; thus, M increases monotonically from 



Off-State, to LS, to Intermediate State, to HS, to VHS (cf. Esin et al. 1997 ). Our results show that, when 
a BHC evolves between LS and HS, it should cross the spectral region in which we observe the VHS (being 
intermediate to the two), and that while doing so, the noise properties of the BHC should change from 
band-limited to power-law. This has been observed recently: Cyg X-1 exhibited 3-12 Hz QPO during its 



1996 transition, when its 2-20 keV intensity was between that of observed HS and LS(Cui 1997); and GRO 



J1655-40, which exhibited 6.5 and 0.8 Hz QPO (Mendez et al. 1998) 



While Cyg X-1 underwent its 1996 transition, QPO (of a broad type, Q^l) were observed, with 
centroid frequencies between 3-12 Hz ( |Cui 1997 ). In Fig. |lj, we show the dependence of Cyg X-1 QPO Vc 



vs. spectral hardness (from Table 2 and 3 of Qui 1997, and W. Cui, priv. comm.). The hardness ratio of 



13.1-60 keV/2.0-6.5 keV was measured by i?XTE/PCA, and cannot be directly compared to the Wide Ratio 
measured by Ginga, although we can reasonably expect that the two are monotonically related. 

Those observations with hardness ratio >0.11 had a type of band-limited noise, with a power-law 
slope '^0.3-0.8 below ~0.5 Hz, then constant between ^0.5-2.0 Hz, breaking again into a power-law of slope 
~ 2.0 above that; during this period (called the "transition phase" by Cui et al. ) the QPO centroid 
frequency was anti-correlated with spectral hardness. During those observations with hardness ratio <0.11 
(called "soft state" by Cui et al. ) the source variability was described as a broken power-law (of slopes ^ 
1 and 2, breaking near i/ ~ 5-15 Hz); in addition, the QPO centroid frequency was below that observed in 
the softest (but relatively spectrally harder) "transition phase" observations, indicating an anti-correlation 
with spectral hardness. 

Comparison between the dependence on the QPO centroid frequency in Cyg X-1 (Fig. with that 
of GX 339-4 and GS 1124-68 (Fig. ||a) shows that in all three objects, the centroid frequency is correlated 
with spectral hardness, which changes to an anti-correlation above some spectral hardness; in addition, all 
three objects show power-law PDS behavior when the centroid frequency is correlated with spectral hardness, 
and band-limited PDS behavior when the centroid frequency is anti-correlated with spectral hardness. 

The differences between the behavior of Cyg X-1 and that which we observe from GX 339-4 and 
GS 1124-68 is that the QPO (in Cyg X-1) are considerably broader QPO (Q~l compared with Q'^10), and 
the 2-20 keV flux is only moderately greater than its Low State flux, and still below that of its High State 



( Pclloni et al. 19971 ), while GX 339-4 and GS 1124-68 were observed with ~ 2-20 keV intensities greater 
than periods where "High State" behavior was observed. It was on the basis of this discrepancy in 2-20 keV 
intensity in Cyg X-1 (pui 1997 ), as well as during a period when QPO were observed at low intensity in 



GS 1124-68 (Belloni et al. 1997), that analogy with the Very High State was discarded. However, the change 
in dependence of the QPO centroid frequency on spectral hardness, associated with the change in the BBN 
variability is suggestive that the phenomena of Cyg X-1 observed during transition between Low and High 
state is the same as that which we describe of GS 1124-68 and GX 339-4 in their Very High State, in spite 
of the differences in 2-20 keV intensity. 
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We conclude that the conditions for transition between the two states occur at luminosities both 
above (in the case of the Very High State) and below (in the case of the Intermediate State) the High-State 
in BHCs. 

We cannot exclude that the assumption that the 2-20 keV flux (usually thought to correspond to 
M) is a monotonic indicator of states is not correct. For example, GS 2023-1-33 exhibited spectral and timing 
behavior which has been classified as "Low State" (spectrally hard, with band- limited noise), even at its 
highest luminosity - although much caution must be employed in interpreting this behavior, as GS 2023+33 



was a peculiar X-ray transient in more than one way (cf. Tanaka fc Lewin 1995 ) 



The strength of the underlying BBN variability is (in GS 1124-68) on average decreased by the 
presence of QPO. One might interpret this to mean that, while the physical process producing the VHS 
QPO is distinct from that of the BBN variability, there is nonetheless interaction between the two processes 
which affects the strength of the BBN variability. 

It does not follow naturally from any other observed behavior that, as a function of spectral hardness, 
the variability ratio should be low (~1), increase with increasing spectral hardness toward the transition 
point (Wide Ratio ~ 0.055, where it is ^J4), and then decrease with increasing spectral hardness. This 
therefore represents a significant and interesting phenomenon, which requires further discussion. It implies 
that the spectrum of the QPO is (in the Low State) similar to the spectrum of the BBN variability; as the 
source spectrum evolves toward the High State, the QPO spectrum becomes spectrally hard, relative to the 
BBN variability, becoming maximally spectrally hard relative to the BBN variability at the same Wide Ratio 
as when the QPO centroid frequency is at its maximum; and then, as the spectrum evolves further toward 
the High State, the spectrum of the QPO again softens toward the spectrum of the BBN variability. (We 
emphasize here that "evolution" is as a function of spectral hardness, not of time). There exists no a priori 
reason that the QPO photon spectrum should at any time be in any way similar to the BBN variability 
photon spectrum - the Variability Ratio could, in principle, take on any value between and infinity. That 
the lowest observed values of the Variability Ratio indicate that the QPO energy spectrum is similar to the 
BBN energy spectrum suggests that the QPO itself is due to the same mechanism as the BBN variability, 
which for some unknown reason becomes enhanced at a particular frequency. 

In summary, the observed dependence on spectral hardness of the BBN variability while QPO 
are observed from GX 339-4 and GS 1124-68, and their position on the color-color diagram spectrally 
intermediate to the Low-State and High-State indicate that the so-called Very High State QPO is due to 
interaction between the physical component responsible for Low State behavior and that responsible for High 
State behavior. Our comparison with QPO observed from Cyg X-1 at lower luminosity, during transition 
between the Low and High States, indicate that these QPO may also be due to such interaction. If this 
inferrence is correct, then the QPO may prove a useful probe in investigating the properties of the intereaction 
of the different accretion components in BHCs. 
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Table 1: Quantity of BHC Data for Timing Analysis, by Source and Temporal Resolution 



Time Res. (sec) 0.00781 0.06250 0.5000 >0.5 Total 



Source 


(ksec) 


(ksec) 


(ksec) 


(ksec) 


(ksec) 


Cyg X-1 


25.2 


4.2 


1.9 


55.8 


87.0 


GS 1124-68 


54.6 


57.9 


268.4 


0.0 


380.8 


GS 2000+25 


30.7 


9.4 


27.1 


71.2 


138.3 


GX 339-4 


15.2 


23.0 


15.2 


18.9 


72.4 



Table 2: QPO Parameters for Best-Fit Cyg X-1 Composite PDS 



Spectral Range 


Energy Range 




FWHM 


%rms 


(Wide Ratio) 


(keV) 


(Hz) 


(Hz) 




Gr 


2.3-4.6 


1.07±0.04 


1.30±0.09 


11.7±0.5 


(0.1-0.15) 


4.6 9.2 


1.15 ±0.03 


1.04 ±0.08 


10.1±0.4 




9.2-18.4 


1.14±0.03 


0.8 ±0.1 


8.2 ±0.6 




2.3-18.4 


1.11 ±0.03 


1.17 ±0.07 


10.8 ±0.4 


G2'' 


2.3-4.6 


(3.14) 


(0.72) 


q r-l-0.3 
o.0_o.4 


(0.15-0.32) 


4.6-9.2 


(3.14) 


(0.72) 


q q-l-0.5 




9.2- 18.4 

2.3- 18.4 


(3.14) 

o 14-^0.09 


(0.72) 
0.7±0.3 


<6.0 

q c-FO.S 
o.0_o.4 



" Fits applied to values in the 0.001 - 64 Hz range 
^ Fits applied to values in the 0.0156 - 64 Hz range 
Values in parenthesis are fixed 
Upper-limits are 3 a 



This preprint was prepared with the AAS lAT^jX macros v4.0. 
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Table 3: %rms of QPO and BBN in GX 339-4 





Energy Range 


QPO %rms 


BBN %rms 


BBN 


FWHM 




Xi (dof) 


Wide Ratio 


(keV) 




(0.001-64 Hz) 


model 


(Hz) 


(Hz) 




0.030-0.043 


2.3-4.6 


1.79±0.04 


3.89±0.09 


PL 


0.53±0.03 


5.80±0.02 


1.36 (139) 


(Fl) 


4.6-9.2 


5.97±0.08 


6.3±0.2 


PL 


0.53±0.02 


5.80±0.01 


1.87 (139) 




9.2-18.4 


10.5±0.2 


5.6±0.5 


PL 


0.55±0.03 


5.78±0.01 


1.36 (140) 




2.3-18.4 


2.89±0.04 


4.45±0.07 


PL 


0.47±0.02 


5.80±0.01 


2.42 (139) 


0.043-0.050 


2.3-4.6 


''•''^-0.08 


13.5±0.8 


DFTPL 


0.59l^„;^ 


6.40±0.04 


2.25 (46) 


(F2) 


4.6-9.2 


8.7±0.1 


23.3±1.5 


DFTPL 


1.07±0.06 


6.37±0.02 


2.60 (46) 




9.2-18.4 


17.4±0.4 


24.3 ±2 


DFTPL 


1.4±0.1 


6.42±0.04 


1.58 (49) 




2.3-18.4 


4.24±0.06 


15.4±0.7 


DFTPL 


0.95±0.05 


6.36±0.02 


4.28 (46) 


0.050-0.060 


2.3-4.6 


1.1±0.2 


8.4±0.2 


FTPL 


(3.58) 


(6.9) 


1.13 (48) 


(F3) 


4.6-9.2 


5.7±0.1 


13.0±0.2 


FTPL 


3.17±0.2 


7.08±0.07 


1.43 (47) 




9.2-18.4 


11.5±0.4 


8.9±0.7 


FTPL 


3.55±0.3 


7.38±0.09 


1.05 (48) 




2.3-18.4 


2.7 ±0.1 


9.5±0.1 


FTPL 


3.40±0.3 


6.9 ±0.1 


1.39 (46) 


0.060-0.075 


2.3-4.6 


1.85±0.06 


8.8±0.1 


FTPL 


(2.27) 


(6.27) 


2.73 (47) 


(F4) 


4.6-9.2 


5.9±0.1 


13.2±0.3 


FTPL 


2.35±0.1 


6.38±0.03 


2.12 (45) 




9.2-18.4 


9.3±0.2 


8.0±0.4 


FTPL 


1.80 ±0.1 


6.39±G.G4 


2.30 (47) 




2.3-18.4 


3.24±0.05 


9.6±0.1 


FTPL 


2.37±0.09 


6.27±0.03 


2.71 (45) 



Fits applied to values in the 0.125 - 64 Hz range 



Table 4: Division of GS 1124-68 Data 



Observation 


Start Time 


End Time 


PDS 






ID 


(UT) 


(UT) 


Freq Range 


<Widc Ratio> 


1 


91/011 00:51 


01:17 


0.0156-8.0 Hz 


0.155 


±0.004 


2 


91/011 19:37 


20:42 


0.0156-8.0 


0.097 


±0.006 


3 


91/011 20:49 


91/012 00:05 


0.0156-64.0 


0.092 


±0.003 


4 (not used) 












5 


91/014 19:48 


23:25 


0.0156-8.0 


0.054 


±0.004 


6 


91/016 18:46 


23:27 


0.0156-8.0 


0.026 


±0.005 


7 


91/017 18:05 


21:57 


0.0156-8.0 


0.014 


±0.001 


8 


91/018 17:30 


21:27 


0.0156-8.0 


0.024 


±0.003 


9 


91/020 15:17 


17:15 


0.0156-8.0 


0.027 


±0.001 


10 


91/022 17:18 


20:09 


0.0156-64.0 


0.058 


±0.004 


11a 


91/025 12:48 


13:10 


0.0156-64.0 


0.036 


±0.001 


lib 


91/025 13:46 


13:56 


0.0156-64.0 


0.044 


±0.001 


12 


91/036 10:34 


11:07 


0.0156-64.0 


0.027 


±0.002 


13 


91/037 10:05 


12:13 


0.0156-64.0 


0.029 


±0.001 


14 


91/044 05:04 


08:38 


0.0156-8.0 


0.024 


±0.004 


15 


91/045 06:11 


08:10 


0.0156-64.0 


0.013 


±0.001 


16 


91/051 23:29 


91/053 04:32 


0.0156-64.0 


0.0093 


±0.0017 


17 


91/056 20:56 


91/058 00:19 


0.0156-64.0 


0.0068 


±0.0007 


18 


91/067 14:28 


91/109 20:13 


0.0156 64.0 


0.0009 


±0.0007 


19 


!)i/i37 03:12 


l)i/20i 21:01 


0.0150 Gl.O 


0.3995 


=0.1578 
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Table 5: GS 1124-68 QPO Fundamental Parameters 



Observation 


Energy Range 




FWHM 


% rms 


xt (dof) 


ID 


(keV) 


(Hz) 


(Hz) 






1 


2.3-4.6 


3.004±0.008 


0.37±0.1 


2.4 ±0.1 


1.61 (39) 




4.6-9.2 


2.98 ±0.01 


0.36±0.05 


4.9 ±0.2 


1.24 (39) 




9.2-18.4 


3.004±0.008 


0.29±0.05 


7.0 ±0.1 


2.46 (39) 




2 3—18 4 


9 QQ_l_n nm 


n i^n-i-n nfi 

U. OUiU.UD 


A n -t-n 1 




2 


2.3—4.6 


no -i-n n 
o.UzitU.Uo 


1 no-i-n 9 


1.6 ±0.2 


1 on i A'\\ 




4.6—9.2 


O.U / ihU.Uz 


n on_Ln n A 
U.80±U.U4 


4.8±0. 1 


1 rtCi / /I Q ^ 
1.99 ^43) 




9.2—18.4 


O.UZihU.Uz 


n £;n_Ln n A 
U.d9±U.U4 


7.9±0.2 


1 rtc\ i A'2\ 
1.99 ^43) 




2 3—18 4 




u.ooiu.uo 


^.oDmu.uo 


1 sn M 1 ^ 


3 


2.3—4.6 


o.ooitU.Uz 


n Q7_i-n riK 
u.y / ilu.ud 


1 7/1 4-n n 
i. i 4itu. Uo 


1 'i'i ( KC\r\\ 
i.oo l^oUUJ 




4.6—9.2 


O.OZOltU.UUo 


u.oyitu.uo 


5.4 ±0.1 


i.o4 I^OUUJ 




9 2—18 4 


^T-l-n nnQ 
o.ooiu.uuy 




1 n fi-t-n 1 

lU.DItU. 1 






2 3—18 4 


o.ozy^u.uu ( 




1 Q-{-n Ci^ 


1 .ou l^OUl ) 











<.U.oU 


1 .00 t^^o J 




4.6—9.2 


D.ozitU.uy 


. 7+0.4 


2.4±0.4 


1.64 (41) 




9 2—18 4 






9 2±0 2 


1 97 (46) 




2.3-18.4 


7.4±0.1 


2.45±0.3 


1.3 ±0.14 


2.29 (41) 


8 


2.3-4.6 


(5.65) 


(2 2) 


0.9±0.07 


1.65 (47) 




4.6-9.2 


(5.65) 


(2.2) 


2.7±0.1 


1.95 (47) 




9.2-18.4 


(5.65) 


(2.2) 


8.8±0.6 


2.00 (47) 




2.3-18.4 


5.65±0.15 


n n + 0.9 
■^■■^ — O 6 


1.0±0.04 


0.99 (42) 


10 


2.3-4.6 


(7.72) 


(3.7) 


0.7±0.1 


1.64 (46) 




4.6-9.2 


8.11 ±0.05 


3.9 ±0.2 


6.0 ±0.3 


1.68 (44) 




9.2-18.4 


8.43 ±0.03 


3.27±0.08 


16.4±0.2 


1.73 (45) 




2.3-18.4 


7.72 ±0.04 


3.7 ±0.1 


2.55 ±0.06 


2.27 (44) 


11a 


2.3-4.6 


(7.5) 


(2.11) 


<0.36 


1.27 (52) 




4.6-9.2 


(7.5) 


(2.11) 


1.4±0.2 


1.75 (50) 




9.2-18.4 


(7.5) 


(2.11) 


<5.1 


0.76 (52) 




2.3-18.4 


7.8±0.3 


2.11±0.6 


0.56±0.07 


2.55 (52) 


lib 


2.3-4.6 


(5.82) 


(0.35) 


1.10±0.08 


1.47 (48) 




4.6-9.2 


5.816±0.009 


0.33±0.08 




1.06 (46) 




9.2-18.4 


5.81±0.01 


0.36 ±0.15 


15.7^1 


0.72 (47) 




2.3-18.4 


5.813±0.008 


0.36±0.07 


3.2±0.3 


1.36 (46) 


15 


2.3-4.6 


4.5 ±0.1 


1 3+"-^ 


0.67±0.05 


1.20 (46) 




4.6-9.2 


4.83±0.09 


1 1 + 0.3 


1.6±0.2 


1.22 (39) 




9.2-18.4 


(4.68) 


(1.14) 


<4.8 


1.58 (51) 




2.3-18.4 


4.6S±0.04 


0.96 ±0.1-}, 


0.75±0.04 


1.36 (39) 



Frequency ranges used to fit and Wide Ratio are in Table 
Values in parenthesis arc fixed 
Upper-limits are 3 tr 



Table 6: QPO Parameters for GS 2000+25 



QPO 



Energy Range 




FWHM 


%rms 


Xi (dof) 


(kcV) 


(Hz) 


(Hz) 






2.3-4.6 




1 0+0.6 
-^■■^-0.3 


4 9+1-^ 


3.40 (29) 
0.79(26) 


4.6-9.2 


(2.63) 


(1.68) 


5.6 ±0.7 


1.52 (28) 


9.2-18.4 


(2.63) 


(1.68) 


5 4+1-^ 


1.91 (30) 


2.3-18.4 


2.63±0.06 


1.7±0.2 


6.3±0.5 


1.39 (26) 



Fit is of single observation period with Wide Ratio— 0.09-0. 1 
Values in parenthesis are fixed 

Fits applied to values in the 0.125 - 8.0 Hz range 
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Fig. 1. — Color-Color Diagram. A composite color-color diagram of all 10 BHC sources observed by 
Gmga (1354-64, 1826-24, 1630-47, LMC X-1, LMC X-3, GS 2000-^25, GS 2023+33, GS 1124-68, Cyg X-1, 
and GX 339-4), which shows the range of spectral hardness observed from each source. Data has been 
corrected for background counts as well as for aspect and deadtime, prior to producing the hardness ratios. 
Each point represents 128sec integration. Only values of the hardness ratio which are >Aa significant are 
shown. 

Fig. 2. — Very High State and Low State QPO, Located on the Wide Ratio vs. Hard Ratio 
CCD. The position on the Wide Ratio vs. Hard Ratio CCD, when QPO is observed from GS 1124-68 (Very 
High State), GX 339-4 (Very High State and Low State from Grebenev et al. 1991), Cyg X-1 (Low State), 
and GS 2000-1-25 (Low State). For clarity, only a few representative points for each source are shown. This 
figure may be directly compared with the full composite CCD in Fig. ^ While the black hole candidates 
occupy four orders of magnitude in Wide Ratio parameter space, the QPO are observed largely within one 
order of magnitude in Wide Ratio. 

Fig. 3. — Source Hardness Ratios. The Wide Ratio of all 10 BHC sources in the present study, which 
shows the range of spectral hardness observed from each source. Each point represents 128 sec integration. 
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Fig. 4. — Intrinsic Source Energy Spectrum vs. Wide Ratio. The intrinsic source photon spectrum 
is indicated by a combination of the ratio of the power-law photon flux to the black body photon flux (2-20 
keV), the photon power-law slope (a=0.9 to 2.9, spaced by 0.4), and the black-body temperature (0.1, 0.25, 
0.4, 0.55, 0.7, 0.85, 1.0, and 1.4 keV). The Wide Ratio is obtained by folding the intrinsic source photon 
spectrum through the Ginga LAC response matrix. 



Fig. 5. — Cyg X-1 Composite PDS, Gl and G2. 8 PDS of Cyg X-1, in 4 energy ranges (indicated in 

each panel). There arc two PDS for each energy range: one from data taken when the source Wide Ratio 
is in the range 0.1-0.15 (Gl; connected by a line to guide the eye), and one from data taken when the 
source Wide Ratio is in the range 0.15-0.32 (G2, unconnected points). The la error bars are included. The 
spectrally soft PDS shows a QPO feature near 1 Hz, while the harder PDS does not have this feature. The 
PDS show roughly equal %rms values above 10 Hz, and %rms values which are not equal below 10 Hz. 
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Fig. 6.— Detailed HID and PDS of GX 339-4 with QPO Each column is headed with the energy 
range of the data used for those PDS, and each row is labeled acicording to spectrum, corresponding to the 
divisions in the HID at the top of the figure. Upper left-hand panel is the Wide Ratio vs. time, during 
the period in which QPO were observed from GX 339-4 (Sept 1988; time=0 is 88/247 13:37 UT). Upper 
right-hand panel is the HID, Intensity Range count-rate vs. Wide Ratio during the Sept 1988 observation. 
The panel has been divided by Wide Ratio value, and each region is labeled. Bottom Panels: The PDS in 
three energy ranges plus the intensity range (labeled by column) of the Sept 1988 observations. Upper-limits 
are 2c7. The data was divided by spectral hardness (each row is labeled, and corresponds to a label in the 
HID above). The BBN PDS characteristics evolve - as the spectrum softens, the BB PDS changes from a 
flat-top-power-law to a pure power-law, while the QPO remains throughout.). 
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Fig. 7.— Detailed HID and PDS of GS 1124-68 with QPO Upper left-hand panel is the Wide 
Ratio vs. time, during the first 40 days of outburst. Upper right-hand panel is the HID, Intensity Range 
count-rate vs. Wide Ratio. The panel has been divided (identicaUy to Fig. |^ by Wide Ratio value, and 
each region is labeled. Bottom Panels: The PDS in three energy ranges plus the intensity range (labeled 
by column) of the Sept 1988 observations. Upper-limits are 2a. The data was divided by spectral hardness 
(each row is labeled, and corresponds to a label in the HID above). The BBN PDS characteristics evolve 
- as the spectrum softens, the BE PDS changes from a flat-top-power-law to a pure power-law, while the 
QPO remains throughout.). 

Fig. 8.— GS 2000+25 PDS of State with QPO. The PDS in the three energy ranges and the sum 
energy range (indicated in each panel). Corrections for background and deadtime have been applied to the 
Power. Upper-limits are 2a. 

Fig. 9. — QPO Centroid Frequency (j/c) and FWHM vs. Wide Ratio. Panel a: Centroid frequency of 
the fundamental QPO, as a function of Wide Ratio. Centroid frequency is maximum near Wide Ratio=0.055. 
Panel b. The FWHM of the fundamental vs. Wide Ratio. The FWHM is greatest near Wide Ratio=0.055. 
Icrerror bars are included on all points. The Cyg X-1 G2-data are not depicted (due to their high Wide 
Ratio) but follow these trends. 
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Fig. 10.^ QPO and BBN %rms, Very High State QPO of GS 1124-68 & GX 339-4. The %rms 
variability in the QPO and BBN components when QPO is present or absent, in GS 1124-68 and GX 339-4. 
Upper-hmits are 3-signia. Points for GS 1124-68 are marked with observation IDs from Table ^; points from 
GX 339-4 are marked with spectral indicator IDs taken from Fig. ^ (a): QPO %rms vs. Energy range from 
GS 1124-68. (h): BBN %rms (when QPO is present) vs. Energy range from GS 1124-68. (c): BBN %rms 
(when QPO is absent) vs. Energy range from GS 1124-68. (d): QPO %rms vs. Energy range from GX 339-4. 
(e): BBN %rms (when QPO is present) vs. Energy range from GX 339-4. 



Fig. 11.— %rms of Low-State QPO vs. Energy Range, Cyg X-1 and GS 2000+25. CygX-1 (open 
circles) and GS 2000+25 (solid squares) Error bars on the %rms are Icr, while on energy they represent the 
energy range of the PDS. The QPO %rms is soft (in Cyg X-1), and flat (in GS 2000+25). This contrasts 
with QPO in GS 1124-68 which, in the same spectral range (Obs. #1-3), as well as in all other obseravations 
and in GX 339-4 as well, is spectrally hard. 



Fig. 12. — Variability Ratio vs. Spectral Hardness . The Variability Ratio of QPO %rms (Eq. |i|) as a 
function of Wide Ratio, from the four QPO sources in the present study. A Variability Ratio of X indicates 
the QPO is X times as strong relative to the BBN variability in the High Energy Range as it is in the Low 
Energy Range. 



Fig. 13. — Variability Ratio vs. QPO Centroid Frequency. The normalized ratio of QPO %rms 
(Eq. |l|) vs. the QPO centroid frequency of sources in the present study. A Variability Ratio of X indicates 
the QPO is X times as strong relative to the BBN variability in the High Energy Range as it is in the Low 
Energy Range. 



Fig. 14. — Cyg X-1 1996 Transition, QPO Centroid EYequency vs. Spectral Hardness. The Vc 

are taken from Table 2 and 3 of Cui et al. (1997), and spectral hardnesses by W. Cui (priv. comm.). The 
value of the spectral hardness ratio is as measured by i?XTE/PCA and cannot be directly compared to the 
Wide Ratio measured by Ginga although the spectral hardness is monotonic, and so can be qualitatively 
compared to the i^c vs. Wide Ratio of GX 339-4 and GS 1124-68 (Fig. ^). Each point is marked with the 
observation number taken from the Table of Cui et al. (1997). 
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Cyg X-1: Ratio = 0.1-0.15 vs. 0.15-0.32 
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